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Stephan’s Quintet (SQ, co-moving radial distance = 85 + 6 Mpc, taken from the NASA/ 
IPAC Extragalactic Database (NED)’) is unique among compact groups of galaxies” ”. 
Observations have previously shown that interactions between multiple members, 
including a high-speed intruder galaxy currently colliding into the intragroup 


medium, have probably generated tidal debris in the form of multiple gaseous and 
stellar filaments®*’, the formation of tidal dwarfs”’*» and intragroup-medium 
starbursts”, as well as widespread intergalactic shocked gas*’°""”. The details and 
timing of the interactions and collisions remain poorly understood because of 

their multiple nature’®”’. Here we report atomic hydrogen (H 1) observations inthe 
vicinity of SQ witha smoothed sensitivity of lo = 4.2 x 10° cm” per channel (velocity 
bin-width Av = 20 kms“; angular resolution = 4’), which are about two orders of 
magnitude deeper than previous observations®***°”. The data show a large HI 
structure (with linear scale of around 0.6 Mpc) encompassing an extended source 

of size approximately 0.4 Mpc associated with the debris field anda curved diffuse 
feature of length around 0.5 Mpc attached to the south edge of the extended source. 
The diffuse feature was probably produced by tidal interactions in early stages of the 
formation of SQ (>1 Gyr ago), although it is not clear how the low-density HI gas 

(Ny,< 10'8 cm”) can survive the ionization by the intergalactic ultraviolet background 
onsucha long time scale. Our observations require a rethinking of properties of gas 
in outer parts of galaxy groups and demand complex modelling of different phases of 
the intragroup medium in simulations of group formation. 


Atomic hydrogen (H 1) is the least bound component of galaxies 
and is therefore the easiest (and hence first) to be stripped off 
and spread around during interactions. Thus, the distribution of 
the very diffuse H Iand its velocity field can provide new information 
about th earliest interactions. To study the diffuse HI associated with 
Stephan’s Quintet (SQ), we carried out deep mapping observations of 
the 21 cm HI emission over a region of around 30’ x 30’ centred onSQ 
(Fig. 1) using the 19-beam receiver of Five-hundred-meter Aperture 
Spherical Telescope (FAST). The FAST observations and data reduc- 
tion are described in the Methods. As shown illustratively in Fig. 2, 
the final data cube includes 304 spectra in Av = 20 kms“ channels 
covering the velocity range of 4,600-7,600 kms“, with an average 
root-mean-square (r.m.s.) noise of 0.16 mJy per beam and an aver- 
age beam size of 2.9’. The mapping satisfies the Nyquist sampling 
criterion with beams separated by 1.4’ in the right ascension direction 
and 1.2’ in the declination direction. The original data cube has aH I 
column-density sensitivity of lo =1.2 x 10” cm” per channel, which 


is improved to lo = 4.2 x 10° cm” per channel when smoothed to 4’. 
Results of the analysis of the whole data cube will be presented 
elsewhere. In this Article we report only the discovery of a large HI 
structure in the velocity range of 6,550-6,750 kms”. 

Figure 3a presents the integrated H I emission map in the veloc- 
ity range of 6,550-6,750 kms overlaid on the deep MegaCam 
optical colour image”. The map has an angular resolution of 
4.0’ and Av=200 kms“ (ten times the channel width) with an 
H 1 column-density error of lo=1.34 x 10” cm”. The base con- 
tour starts from N,,,= 7.4 x 10” cm”? (at a 5.50 level). The map 
shows a large H I structure of around 0.6 Mpc in size, which has 
two parts: an extended source centred on SQ and a diffuse fea- 
ture attached to the south edge of the source. The extended 
source encompasses the previously detected 6,600 kms?HI1 
component associated with the debris field®’. As shown by the cyan 
contours in Fig. 3a, the high-resolution (beam = 19.4” x 18.6”) Very 
Large Array (VLA) observations® detected only the high- density 
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Fig. 1| Member galaxies plus main tidal features in SQ and the sky coverage 
of the FAST observations of the SQ field. The background image is the 
inverted greyscale map of the deep Canada-France-Hawaii Telescope (CFHT) 
MegaCam r-band image”. The NGC name, radial velocity and the Hubble type 
are provided for each member galaxy inSQ. The core members include NGC 
7317, NGC 7318a and NGC 7319, which are located near the group centre and 
have similar radial velocities (6,680 + 85 kms‘). NGC 7320cis an ‘old intruder’, 
which may be responsible for the formation of the inner and outer tails®”. 
NGC 7318b is a ‘new intruder’, whichis currently colliding into the intragroup 


part of the 6,600 kms? component (N,,=5.8 x 10” cm”), which is 
confined to the central region (D = 0.1 Mpc) of the extended source. 
Most of the high-density H I gas traces the optically detected inner 
and outer tails® plus a compact cloud (north-west of NGC 7319) 
coincident with the intragroup-medium starburst SQ-A”. Single-dish 
H I mapping observations by the Arecibo Telescope and the 
Green Bank Telescope, which detected lower density H I gas at 
Nu, = 5 x 10'8 cm” albeit with lower angular resolutions (>3’), have 
found evidence for this component to be extended ona scale of 
around 0.2 Mpc (refs. *°”"). Our deeper FAST map shows an even 
larger size with a diameter of around 0.4 Mpc. The diffuse feature 
has a characteristic column density of approximately 7 x 10” cm? 
in an elongated and curved structure of around 0.5 Mpc in length. 
It appears at the bottom of the FAST mapping and therefore may 
well reach beyond the map. The faint optical halo (the yellowish diffuse 
light around SQ in the optical colour image) discovered previously” 
lies inside the extended source and has no spatial overlap with 
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medium of SQ and triggering a large scale shock®. NGC 7320 isa foreground 
galaxy. The small blue circles mark the positions of individual beams inthe HI 
mapping observations by the FAST 19-beam receiver. The observations were 
carried out at 16 slightly separated pointingsina4 x 4 rectangular grid. The 
small red circles mark the central positions of the FAST 19-beam receiver in 
these pointings, and the characters P1, P2, ..., P16 inside the circles identify the 
different pointings. The large magenta circles (D = 2.9’, that is, the half-power 
beam size) show the coverage of the 19 beams in the first pointing (see Methods 
for more details). 


the newly discovered diffuse feature. The first moment map in 
Fig. 3b shows that in the velocity field the diffuse feature is linked 
smoothly with the extended source. The two green boxes marked 
by characters A and B in Fig. 3a cover the entire diffuse H I feature. 
The sum of all spectra in these two boxes provides a good meas- 
ure of the spectrum of the diffuse feature, which is presented in 
Fig. 3c. The spectrum has a flux-density-weighted mean velocity of 
6,633 kms‘and a rather narrow line width of AV) = 160 km s 7 (AV39 
is the line width measured at 20% of the peak). The integrated 
flux is 0.42 + 0.03 Jy kms, corresponding to an H I mass of 
(7.14 0.5) x 10°Mo, which is only around 3% of the total HI mass of 
SQ (2.45 x 10°°M)”. Itis worth noting that, although the Green Bank 
Telescope mapping observations found 65% more HI than the VLA 
observations, SQ is still slightly deficient in H 1 abundance compared 
with normal galaxies (by a factor of around 1.3)*”. The very diffuse 
H1(Ny,<3 x 10'8 cm”) discovered in this work does not change this 
HI deficiency significantly. 


Fig. 2| Illustrative plot of the 304 spectra of the H1 emission in the SQ field. 
The underlying optical colour image (u, g, r) is obtained in the deep CFHT 
MegaCam observation with limiting surface brightness of 29.0, 28.6 and 

27.6 mag arcsec ’ for the three bands, respectively”. The spectra cover the 
velocity range of 4,600-7,600 kms‘ withan average r.m.s. of 0.16 mJy per 


Two new detections of unresolved sources can also be foundin Fig. 3a. 
NGC 7320a, detected with the signal-to-noise ratio S/N = 36, hasan HI 
mass of M,,, = (6.3 + 0.2) x 10°M, anda vy,, = 6,702 + 24 kms‘. The other 
source Anon7,a 4.40 detection, has an HI velocity of 6,654 +16kms* 
and an HI mass of (2.2 + 0.5) x 10°M,. More discussions about these 
two sources are given in the Methods. 

We examine in Fig. 4 the individual spectra in box A and box B to 
investigate the physical nature of the diffuse feature. No diffuse stel- 
lar radiation is detected in these regions down to the limit of the deep 
MegaCam image. Spectra of beams with detections of S/N > 4 are 
marked by pink boxes and those with 3 < S/N < 4 by green boxes. Gal- 
axies brighter than the r-band magnituder = 20 mag found inthe SDSS 
photometric redshift (photo-z) catalogue” are also marked in Fig. 4. 
Only two of them have photo-z < 0.1 (marked by red circles) and the 
remaining 28 have photo-z = 0.1 (orange circles). Given the loerror of 
the photo-z (6z/(1+ z) = 0.02)” and the redshift of SQ (z= 0.02), gal- 
axies with photo-z = 0.1 are very unlikely to be at the same redshift 
of SQ (probability < 0.001). We can rule out with high confidence the 
possibility of the diffuse feature being associated with a collection of 
gas-rich galaxies (evenincluding those fainter thanr = 20 mag), because 
it needs at least four such galaxies to cover all beams with significant 
detections of S/N > 4 (one for those in the upper-left corner of box A, 


beam per channel (Av=20 kms‘). The centre of each spectrum coincides with 
the pointing position of the FAST beam with which it was obtained. The final 
data cube includes 304 spectrain20 kms ‘bins. The red circle inthe 
bottom-left corner shows the size of the FAST beam (2.9, or 72.5 kpc onalinear 
scale). 


two for those in the right half of box A and one for that in box B) and 
the probability that they happen to have about the same radial velocity 
is extremely low. 

H I clouds without stellar counterparts have been found in and 
around many galaxy groups/clusters”*°. Most of them, as the authors 
of ref. 7” have argued, can be explained by tidal debris of galaxy inter- 
action involving very extended HI disks instead of ‘dark’ or ‘almost 
dark’ galaxies. Given its location and velocity, the diffuse feature is 
most likely to be related to the debris field. A hypothetical scenario for 
the formation of the diffuse feature is that NGC 7320a (v= 6,702 kms 
and currently around 300 kpc away from the SQ centre) passed 
through the SQ centre approximately 1.5 Gyr ago (assuming a rela- 
tive transverse velocity of 200 kms“) and pulled out from one of 
the core member galaxies of SQ a tidal tail, which developed into 
the diffuse feature we see now. Another possibility is that, like the 
large Leo Ring (D = 0.25 Mpc)’, the diffuse feature could be the prod- 
uct of a high-speed head-on collision between another old intruder 
and one of the core members of SQ. In this scenario, the collision 
triggers an expanding density wave that pushes gas in an extended 
HI disk of the target galaxy outwards to form a very large ring, of 
which the diffuse feature is the high-density part. A candidate for 
such an intruder could be Anon 4 (v= 6,057 kms7, M,,,= 1.1 10°M,)8, 
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Fig. 3 | The H1emission in the velocity range of 6,550-6,750 km s*. 

a, Contour map of the integrated HI emission in the velocity range of 
6,550-6,750 kms ‘ overlaid on the composite colour image (u, g, r) of the 

deep CFHT MegaCam observation. The red circle at the bottom right illustrates 
the angular resolution (half-power beam size = 4.0’) of the FAST map after 
smoothing. The contours start from Ny, = 7.4 x 10” cm? (at 5.5olevel) with an 
increment ofa factor of 2. The red lines delineate the boundary of the FAST 
observations. The cyan contours inthe centre are adopted fromthe VLA 
observations for the 6,600 kms‘ component of SQ$, with angular resolution 

of 19.4” x 18.6”. They have the base level at N,,, =5.8 x 10” cm and the increment 
of a factor of 2. The area occupied by the newly discovered diffuse feature is 


which spatially coincides with optical galaxy LEDA 141041 (B band 
magnitude B = 18.4 mag). It has a relative radial velocity of around 
600 kms ‘anda projected distance of approximately 0.2 Mpc from 
the SQ centre. If the relative transverse velocity is around 200 kms", 
it would have taken approximately 1 Gyr for Anon 4 to move to the 
current position after the collision. Both scenarios proposed above 
suggest a formation time of the diffuse feature of more than 1 Gyr ago. 
They are both based on analogies to cases studied in simulations in 
the literature, which demonstrate that diffuse H1 features without a 
stellar component can be produced in galaxy-galaxy interations”””*. 
However, two questions remain to be answered: (1) Can the tidal fea- 
ture in either of the scenarios survive the subsequent interactions 
that triggered the formation of the inner and outer tails of SQ about 
(3-8) x 108 years ago'®!°? (2) Can H I structures with column density 
as low as N,,, < 10'° cm” exist on timescales of around 1 Gyr? These 
questions can only be answered by more sophisticated models that 
are built upon the existing simulations for the formation and evolution 
of SQ’. It has been argued that cold gas of N,,, <2 x 10” cm™® cannot 
stay neutral in the intergalactic ultraviolet background radiation for 
more than 500 Myr (refs. "”°). A plausible solution for this problem is 
the physical mechanism involving the transition between ionized and 
neutral phases due to thermal instabilities in the low-density gas*°”). 
New simulations, which are beyond the scope of this Article, shall 
explore this mechanism. 
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filled with a hatch pattern consisting of thin yellow dotted lines. The two 
green boxes marked by characters A and B cover the diffuse feature. b, False 
colour map of the velocity field of the HI emission in the velocity range of 
6,550-6,750 kms ‘overlaid by the same contour map shownina. The units 
of thexand yaxes are in pixels (pixel size: 1.4’ x 1.2’), and the colour of every 
pixel represents the flux-density weighted mean velocity whose scale is 
giveninthescale bar onthe right. c, The HI spectrum of the diffuse feature. 
The red line is the summation of all spectra inside box A and the blue line the 
summation of those inside box B. The black line is the sum of the two spectra 
of boxes A andB. 
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Methods 


FAST observations 

The deep HI mapping observations were carried out in September and 
October of 2021 using the FAST 19-beam receiver in the standard ON- 
OFF mode with the total observation time of 22.4 hincluding overheads 
(Extended Data Table 1). The FAST 19-beam L-band Array is currently the 
largest multibeam feed array for HI observations in the world. Details 
about its properties and performance can be found in ref. ”. The 19 
beams are arranged ina hexagonal configuration with the neighbour- 
ing beams separated by 5.7’. The observations have a central frequency 
of 1,391.64 MHz and a frequency coverage of 1,050-1,450 MHz witha 
resolution of 7.63 kHz (Av =1.65 kms“). For the 19 beams, the average 
half-power beam-width at 1,391 MHzis 2.9’ (Extended Data Table 2). To 
meet the Nyquist sampling criterion and fill the gaps between beams in 
the focal plane, we carried out 16 pointings ina 4 x 4 rectangular grid 
inthe north-up orientation (Fig. 1). The final mapping covers a region 
of around 30’ x 30’ centred on SQ with 304 sky pixels (beam positions 
inthe sky), and the separation between the nearest pixels is 1.4’ in the 
right ascension direction and 1.2’ in the declination direction. The 1o 
pointing error of individual beams is 7.9” (ref. ”). At each pointing, six 
cycles of ON-OFF integrations were conducted, with the OFF position 
at 40’ southeast from the ON position. Each ON or OFF took 300 sinte- 
gration with asampling frequency of 1 Hz. The total on-target time for 
each pixel was 1,800 s (Extended Data Table 1). To minimize the effects 
of standing waves and sidelobes, all observations were confined within 
the zenith angle of 20°. 

Compared with the scan-mapping mode, which is more suitable 
for large sky surveys such as the FAST Extragalactic HI Survey”, our 
observational strategy provides an alternative for deep and small maps 
(<30’), which can take advantage of the ON-OFF mode in more accu- 
rately removing various systematic effects such as the standing waves 
and baseline wobbling. 


Data reduction and calibration 

For each sky pixel observed by a given beam, we reduced the spectral 
data following a similar procedure to that presented in ref. **. The spec- 
tra of the two polarizations were reduced separately (and eventually 
combined after the consistency check). The data were grouped in ON- 
OFF cycles. Each ON (or OFF) has 300 samplings which were averaged 
and calibrated (that is, converted from digital counts to kelvin), result- 
ing ina single raw spectrum. During the observations, a calibration 
signal of 10 K was injected for a duration of 20 s at the beginning of 
every ON-OFF cycle, and these data were used to calibrate the antenna 
temperature T,. The calibration error is on the order of 10%. Repeating 
this, we obtained a raw spectrum for every ON or OFF. As examples, the 
upper panel of Extended Data Fig. 1 presents the individual raw ON-OFF 
spectra obtained by the MOI beam during the first pointing observa- 
tion (the P1 pixel in Fig. 1). The mean of these spectra is presented in 
the middle panel of Extended Data Fig. 1. It is affected significantly by 
the standing waves, which can be well fitted by a sine function locally. 
The spectrum has been converted from 7, (in K) to flux density (in mJy). 
The gain factor that converts 7, to flux density (in the units of K Jy”) 
depends on the frequency and varies from beam to beam. The values 
for individual beams at 1,391 MHz, derived by interpolating values 
at other frequencies adopted from ref. ”, are presented in Extended 
Data Table 2. The next step was to remove the standing waves together 
with the baseline from the spectrum. The bottom panel of Extended 
Data Fig. 1 presents the final spectrum at sky pixel Plafter the subtrac- 
tion of a baseline modelled by a sinusoidal (representing the standing 
waves) plus a polynomial (for the baseline gradient). We converted the 
frequency to velocity by adopting the optical redshift convention and 
the local standard of rest reference frame and rebinned the spectrum 
into bins of Av=20 kms”. The above process was carried out repeat- 
edly for every sky pixel observed in our observations. 


Itis worth noting that, if the standing waves and baseline vary signifi- 
cantly during a given observation, it is better to carry out the removal 
of the standing waves and baseline for spectra of individual ON-OFFs 
instead of doing it for their mean. However, in a test in which we did 
the standing waves and baseline subtraction for each ON-OFF and 
then used the median of the baseline-subtracted spectra of ON-OFFs 
as the final spectrum at a given pixel, we got anoisier product. It seems 
that the wavelength of the standing waves (in the frequency domain) 
inagiven observation is rather constant (although the phase changes 
fromcycleto cycle). Consequently, the effect of the standing waves in 
the mean of the spectra of individual ON-OFFs is still a well-defined 
sinusoidal with the same wavelength, which can be easily removed. 
Hence, because the mean spectrumis less noisy than individual spectra 
and therefore a more accurate model for the standing waves and the 
baseline can be obtained, subtracting the baseline model from the 
mean spectrum can achieve a better result. 

The final data cube was constructed from the 304 individual 
spectra obtained using this technique, with a velocity coverage of 
4,600-7,600 km s‘in20 kms ‘bins. A uniform half-power beam-width 
of 2.9’ was adopted for all spectra, neglecting the small variation of the 
beam size among different beams. 

The actual half-power beam-width of individual beams at 1,391 MHz, 
derived from interpolations of values adopted from ref. ’, are listed in 
Extended Data Table 2. The units of the flux density of the spectra in 
the cube are mJy per beam. To find the flux density in units of mJy for 
a given spectrum, a factor of A = (B/2.9’)* should be multiplied to the 
value taken from the cube, where B is the half-power beam-width (in 
arcmin) of the beam with which the spectrum was observed. 

Inthe 16-pointing observations, each of the 19 beams covered 16 adja- 
cent sky pixels. The mean and the standard deviation of the measured 
r.m.s. noise of these 16 spectra are also listed in Extended Data Table 2. 
The r.m.s. noise of each spectrum was measured in the two velocity 
intervals of 4,700-5,000 kms tand 7,000-7,500 kms, where no HI 
signal was detected. Beam M16 stands out as the noisiest beam in the 
array, witha meanr.m.s. of 0.26 mJy per beam and a standard deviation 
of 0.14 mJy beam. The false colour map of the r.m.s. noise in the left 
panel of Extended Data Fig. 2 shows that indeed pixels inthe top-right 
corner that were covered by beam M16 have higher noise than others. 
The histogram of the distribution of the r.m.s. noise at all sky pixels 
is shown in the right panel of Extended Data Fig. 2. The mean of the 
r.m.s. is 0.16 mJy per beam with a standard deviation of 0.05 mJy per 
beam. It is worth noting that strong radio frequency interference (RFI) 
in the frequency range of our observations was a serious issue in the 
early stages of this project. The operation team of FAST did excellent 
work in discovering and removing the source of the RFl ina relatively 
short time. All of our observations were carried out after the removal 
of the RFI source. Consequently, our observations were not affected 
by any significant RFI. 


Smoothing 

The H1 data cube obtained above is highly redundant in the sense 
that a sky area of the size of a single beam (D = 2.9') is covered by 
multiple beams (beam separation 1.4' x 1.2'). When making channel 
maps and integrated emission maps from the data cube, applying a 
Gaussian-kernel convolution (that is, smoothing) makes good use of 
this redundancy. This minimizes the noise due to the signal fluctuations 
in adjacent beams and results in significant improvement inthe HI 
column-density sensitivity. The only disadvantage is a slight degrada- 
tionin the angular resolution. For single-channel maps, the meanr.m.s. 
of 0.16 mJy per beam corresponds to aH I column-density sensitivity 
of lg =1.2 x 10” cm” per channel (Av = 20 kms"). When a smoothing 
with a Gaussian kernel of full-width at half-maximum (FWHM) = 2.8’ is 
applied, the HI column-density sensitivity is improved by a factor of 
2.9 to lo= 4.2 x 10° cm” per channel whereas the angular resolution 
is degraded only slightly (by a factor 1.4) to 4.0’. The improvement in 
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the HI column-density sensitivity is particularly important for the 
exploration of diffuse extended emission. In Extended Data Fig. 3 we 
present the contour map of the integrated HI emission in the velocity 
range of 6,550-6,750 kms ‘before the smoothing. Compared with the 
map after the smoothing (Fig. 3a), the low HI column-density features 
in Extended Data Fig. 3 are more fragmented mainly because of the 
signal fluctuations in adjacent beams. Ona linear scale the pre- and 
after-smoothing resolutions are 72.5 kpc and 100 kpc, respectively, 
which is not a significant difference given that we are searching for 
extended diffuse HI gas ona linear scale of a few 100 kpc. 


Sidelobe correction 

The data cube is corrected for the sidelobes using the images of individ- 
ual beams of the 19-beam receiver”, which provide information about 
the point spread functions of the beams. For each beam, a ‘sidelobe 
responsivity function’ is defined by the difference between the point 
spread function and the ‘main beam’, the latter being approximated by 
a two-dimensional Gaussian with the FWHM equal to the half-power 
beam-width. In the calculation of the sidelobe corrections, we con- 
sider only the effects due to the central extended source associated 
with the SQ group. The FAST observations also detected numerous 
other HI sources in the SQ neighbourhood in the velocity range of 
5,500-7,000 kms”. They are much fainter than the central SQ source 
and therefore their contributions to the sidelobes are neglected. 

The first step is to estimate the sidelobe contribution to the map 
of integrated H I emission in the velocity range of 6,550-6,750 kms1, 
which encompasses the peak of the HI spectrum of SQ”. The original 
integrated H I emission map (the observed map) is first deconvolved 
with the main beam of MO1 (the central beam of the 19-beam receiver). 
Then all pixels outside a circular aperture of D = 10.6’, within which the 
central extended source is located, are masked. The result is then taken 
as our approximation for the ‘truth map’. The sidelobe contribution 
to any given pixel in the observed map is estimated by the following 
equation: Fidetove Xi Vj) = Limon! Xin) * Ry Xm Xv In — Yj), Where x,and 
y;are the coordinates of the pixel centre, 7(X,,y,) is the flux in the truth 
map in the pixel at x,, andy, R, is the sidelobe responsivity function of 
the beam that is pointed at the pixel (x, y,), and ),,, and ),, are summa- 
tions along the xandy directions, respectively. Extended Data Figure 4 
presents the map of the sidelobe contribution estimated using this 
method overlaid by contours of the map of integrated HI emission in 
the velocity range of 6,550-6,750 km s ‘(without the sidelobe correc- 
tion, smoothed by a Gaussian kernel of FWHM = 2.8’). It shows that the 
sidelobes contribute significantly at the edge of the debris field but 
have a minimal effect on the diffuse feature in the south. 

Neglecting the frequency dependencies of the shapes of both the 
sidelobes and the central SQ source, we estimate the sidelobe contribu- 
tion to each channel in the data cube by scaling the map in Extended 
Data Fig. 4 witha factor of C, = S,,/S¢,550-67s0, Where S, (in mJy) is the flux 
density of SQin the given channel and S¢ 550-6750 (in mJy km s“) is the inte- 
grated flux of SQin the velocity range of 6,550-6,750 km s”. Finally, we 
obtain the sidelobe corrected data cube by subtracting the estimated 
sidelobe contribution from every channel map in the cube. The 304 
spectrain the resulting data cube are presented illustratively in Fig. 2. 


Detections of two new H I sources in the SQ neighbourhood 

In the velocity range of 6,550-6,750 km swe detected two new unre- 
solved H Isources in the SQ neighbourhood (Extended Data Table 3). 
Their HI spectra are presented in Extended Data Fig. 5. The spectrum 


of the source associated with NGC 7320a shows a typical double-horn 
profile consistent with the highly edge-on optical morphology of the 
galaxy. To confirm the association of the HI source and the optical 
galaxy, we made along-slit optical spectroscopic observation (1 hexpo- 
sure) for NGC 7320a on the night of 21 December 2021 using the 2.4 m 
telescope at Lijiang Observatory. The optical spectrum is presented 
in Extended Data Fig. 6. The radial velocity obtained from the optical 
spectrum is 6,729 + 59 kms“, which is consistent with the H I velocity 
(6,702 +24 kms‘). The other source is a 4.40 detection without obvi- 
ous optical counterpart and therefore we name it Anon 7 following the 
convention in the literature®”, 


Data availability 


Observational data are available from the FAST archive (http://fast.bao. 
ac.cn) lyear after data collection, following FAST data policy. The data 
that support the findings of this study are openly available in Science 
Data Bank at https://www.scidb.cn/s/jilfee. 


Code availability 


The Python and IDL code for the data reduction pipeline is available at 
https://www.scidb.cn/s/jilfee. 
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Extended Data Fig. 1| Examples ofintermediate products of spectral data 
reduction. Upper Panel: Individual ON-OFF spectra obtained by beam MO1in 
Pointing 1 observation (i.e. the data obtained at the pixel Plin Figure 1). An 
offset of 0.2 Kis added to every spectrum relative to the previous one in order 
to make them separated from each other. Middle Panel: The mean of the 
ON-OFF spectra inthe upper panel. The units of the flux density are converted 
from K to mJy using the gain factor taken from Extended Data Table 2. The red 


6000 
Visr [km/s] 


—— Rebin to 20 km/s 


line shows the baseline model derived by fitting the parts of the spectrum free 
of signals (marked by orange color) using a sinusoidal (representing standing 
waves) plus a polynomial (for baseline gradient). Bottom Panel: The spectrum 
after the baseline removal. The frequency is converted to radial velocity inthe 
optical convention relative to the local standard rest reference frame (LSR). 
The orange line presents the rebinned spectrum with the velocity bin-width of 
20kms". 
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Extended Data Fig. 2 | Spectral R.M.S. noise of individual sky pixels. Left: The false color image of the r.m.s. noise at different sky pixels. Right: Histogram of 
the distribution of the r.m.s. noise. 


Extended Data Fig. 3| The Hl emission in the velocity range of 6550-6750 
kms? (unsmoothed). Contour map of integrated HI emission (unsmoothed) in 
the velocity range of 6550-6750 kms‘ overlaid on the composite color image 
(u, g,r) of the deep CFHT MegaCam observation. The red circle at bottom-right 
illustrates the angular resolution (FWHM =2.9') of the FAST map (unsmoothed). 
The contours start from N,,, = 7.4 x 10” cm? (at 2-o level) with an increment ofa 


factor of 2. The red lines delineate the boundary of the FAST observations. The 
cyancontours inthe center are adopted from the VLA observations for the 
6600 component of SQ, with angular resolution of 19.4"18.6". They have the 
base level at N,,, =5.8 x10” cm and the increment of a factor of 2. The same Box 
Aand Box B that mark the location of the diffuse feature in Figure 3a are plotted 
here to facilitate the comparison between the two figures. 
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Extended Data Fig. 4 | Sidelobe contribution to the integrated HI emission the convolution by a Gaussian kernel of FWHM = 2.8'). The contours start from 
in the velocity range of 6550-6750 kms‘. False color image of the sidelobe 20 mJy kms ‘beam (corresponding to Ny = 7.4 x10” cm for the original beam 
contribution overlaid by the contour map of integrated HI emissioninthe of FWHM =2.9') with anincrement of a factor of 2. 

velocity range of 6550-6750 kms ‘(before the sidelobe correction and after 
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Extended Data Fig. 5| HI spectra of NGC 7320a and Anon 7. HI spectra of two newly detected HI sources (NGC 7320a and Anon 7) inthe SQ neighborhood. 
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Extended Data Fig. 6 | Optical spectrum of NGC 7320a. Optical spectrum of NGC 7320a obtained at Lijiang 2.4 meter telescope. The vertical dotted lines mark 
the positions of corresponding emission/absorption lines redshifted according to the best-fit redshift z = 0.02243 (v=6729 kms”). 


Extended Data Table 1| FAST Observations 


Obs. date ttotal R.A.> Decl.» cent freq. ON-OFF ton 

(secs) (2000J) (2000J) (MHz) cycles (secs) 
Pointing 1 4 2021-09-09 5580 22:35:43.37 +33:59:20.6 1391.64 6 1800 
Pointing 2 4 2021-09-14 5580 22:35:50.39 +33:59:20.6 1391.64 6 1800 
Pointing 3 4 2021-09-11 5580 22:35:57.31 +33:59:20.6 1391.64 6 1800 
Pointing 44 2021-09-13 5580 22:36:04.23 +33:59:20.6 1391.64 6 1800 
Pointing 5 2021-09-15 4860 22:35:43.37 +33:58:06.1 1391.64 6 1800 
Pointing 6 2021-09-18 4860 22:35:50.39 +33:58:06.1 1391.64 6 1800 
Pointing 7 2021-09-20 4860 22:35:57.31 +33:58:06.1 1391.64 6 1800 
Pointing 8 2021-09-24 4860 22:36:04.23 +33:58:06.1 1391.64 6 1800 
Pointing 9 2021-09-26 4860 22:35:43.37 +33:56:51.5 1391.64 6 1800 
Pointing 10 2021-09-27 4860 22:35:50.39 +33:56:51.5 1391.64 6 1800 
Pointing 11 2021-09-29 4860 22:35:57.31 +33:56:51.5 1391.64 6 1800 
Pointing 12 2021-09-28 4860 22:36:04.23 +33:56:51.5 1391.64 6 1800 
Pointing 13 2021-09-30 4860 22:35:43.37 +33:55:36.9 1391.64 6 1800 
Pointing 14 2021-10-03 4860 22:35:50.39 +33:55:36.9 1391.64 6 1800 
Pointing 15 2021-10-01 4860 22:35:57.31 +33:55:36.9 1391.64 6 1800 
Pointing 16 2021-10-05 4860 22:36:04.23 +33:55:36.9 1391.64 6 1800 


a: Total observation time including all overheads. 

b: Coordinates of the pointing position of the central beam M01. 

c: Total on-target time. 

d: These earlier observations have higher overheads than the later observations. 
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Extended Data Table 2 | Properties of beams in the 19-beam receiver 


Beam half-power size 4 gain factor > rms © dorms 4 
ID [arcmin] (Jy/K) (mJy beam) (mJy beam) 
M01 2.8335 16.0128 0.18 0.07 
Mo2 2.8635 15.0520 0.14 0.03 
M03 2.8535 15.5324 0.13 0.02 
M04 2.8780 15.2122 0.16 0.04 
M05 2.8480 15.6638 0.15 0.04 
M06 2.9680 15.0520 0.16 0.04 
M07 2.8690 14.9207 0.16 0.04 
Mo8 2.9090 14.3089 0.17 0.03 
Mog 2.8935 14.3089 0.16 0.04 
M10 2.9590 15.2409 0.14 0.03 
M11 2.9290 14.6004 0.14 0.03 
M12 2.9590 14.2802 0.15 0.03 
M13 2.9190 15.2409 0.14 0.03 
M14 2.9635 14.6004 0.14 0.03 
M15 2.8990 13.9599 0.16 0.02 
M16 2.9790 13.3194 0.26 0.14 
M17 2.9290 13.8285 0.17 0.05 
M18 2.9545 13.7710 0.17 0.04 


M19 2.9190 14.2802 0.15 0.03 


a: Half-power beam size at 1391MHz. 

b: Gain factor at 1391MHz. 

ce: Mean r.m.s. of the 16 spectra (bin-width of 20km s") obtained with the beam. 
d: Standard deviation of the r.m.s. 


Extended Data Table 3 | Two new HI sources in the SQ neighborhood 
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Source ID R.A. Decl. VHI Av20 Mui 
(2000J) (2000J) (km s*) (km s*) (108 Mo) 
Sa ee ee ee ee 
NGC 7320a 22:36:32.2 +33:47:46 6702 + 24 360 6.3 +0.2 


Anon 7 22:35:19.3 +34:08:00 6654 + 16 160 2.2+0.5 


eee 


